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Abstract. A gauge singlet scalar with non-minimal coupling to gravity can drive inflation and 
later freeze out to become cold dark matter. We explore this idea by revisiting inflation in the 
singlet direction (S-inflation) and Higgs Portal Dark Matter in light of the Higgs discovery, 
limits from LUX and observations by Planck. We show that large regions of parameter space 
remain viable, so that successful inflation is possible and the dark matter relic abundance 
can be reproduced. Moreover, the scalar singlet can stabilise the electroweak vacuum and 
at the same time overcome the problem of unitarity-violation during inflation encountered 
by Higgs Inflation, provided the singlet is a real scalar. The 2 -it Planck upper bound on Ug 
imposes that the singlet mass is below 2 TeV, so that almost the entire allowed parameter 
range can be probed by XENONIT. 
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1 Introduction 

There is strong evidence in support of the idea that the Universe underwent a period of 
primordial inflation. In particular, the observation of adiabatic density perturbations with a 
spectral index which deviates from unity by a few percent [1] is consistent with the generic 
prediction of scalar field inflation models. However, the identity of the scalar field responsible 
for inflation remains unknown. Another unsolved problem of similar importance for cosmol¬ 
ogy is the nature of dark matter (DM). While it is possible to explain DM by the addition 
of a new particle, there is presently no experimental evidence for its existence or its identity. 

The most studied proposal is that DM is a thermal relic weakly-interacting massive 
particle (WIMP). WIMPs typically have annihilation cross-sections comparable to the value 
required to reproduce the observed density of DM, the so-called “WIMP miracle”. Neverthe¬ 
less, the non-observation of any new weak-scale particles at the LHC beyond the Standard 
Model (SM) places strong constraints on many models for WIMPs, such as in supersymmetric 
extensions of the SM. The absence of new particles may indeed indicate that any extension 
of the SM to include WIMP DM should be rather minimal. In the present work we therefore 
focus on a particularly simple extension of the SM, namely an additional gauge singlet scalar, 
which is arguably one of the most minimal models of DM [2-4]. 
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A similar issue arises from recent constraints on inflation. In fact, the non-observation 
of non-Gaussianity by Planck [5] snggests that the inflation model should also be minimal, in 
the sense of being dne to a single scalar field. The absence of evidence for new physics then 
raises the question of whether the inflaton scalar can be part of the SM or a minimal extension 
of the SM. The former possibility is realized by Higgs Inflation [6], which is a version of the 
non-minimally coupled scalar field inflation model of Salopek, Bond and Bardeen (SBB) [7] 
with the scalar field identified with the Higgs boson. A good example for the latter option 
are gauge singlet scalar extensions of the SM, because the DM particle can also provide a 
well-motivated candidate for the scalar of the SBB model. In other words, in these models 
the same scalar particle drives inflation and later freezes out to become cold DM. 

The resulting gange singlet inflation model was hrst considered in [8] , where it was called 
S-inflation (see also [9]).^ All non-minimally conpled scalar field inflation models based on 
the SBB model are identical at the classical level but differ once quantnm corrections to the 
inflaton potential are included. These result in characteristic deviations of the spectral index 
from its classical value, which have been extensively studied in both Higgs Inflation [6, 11-15] 
and S-inflation [16]. 

Since the original stndies were performed, the mass of the Higgs boson [17] and the 
Planck resnlts for the inflation observables [1] have become known. In addition, direct DM 
detection experiments, such as LUX [18], have imposed stronger bounds on gauge singlet 
scalar DM [19-25]. This new data has important implications for these models, in particular 
for S-inflation, which can be tested in Higgs physics and DM searches. The main objective 
of the present paper is to compare the S-inflation model with the latest results from CMB 
observations and direct DM detection experiments. 

We will demonstrate that — in spite of its simplicity — the model still has a large viable 
parameter space, where the predictions for inflation are consistent with all cnrrent constraints 
and the observed DM relic abnndance can be reprodnced. In addition, we observe that this 
model can solve the potential problem that the electroweak vacuum may be metastable, 
because the singlet gives a positive contribution to the rnnning of the quartic Higgs coupling. 
Intriguingly, the relevant parameter range can be almost completely tested by XENONIT. 

Another important aspect of our study is perturbative unitarity-violation, which may 
be a significant problem for Higgs Inflation. Since Higgs boson scattering via graviton ex¬ 
change violates unitarity at high energies [26, 27], one might be worried that the theory is 
either incomplete or that perturbation theory breaks down so that unitarity is only conserved 
non-perturbatively [28-31]. In both cases there can be important modification of the infla¬ 
ton potential due to new physics or strong-conpling effects. Indeed, in conventional Higgs 
Inflation, the nnitarity-violation scale is of the same magnitude as the Higgs field during 
inflation [14, 32], placing in doubt the predictions of the model or even its viability. 

In contrast, we will show that S-inflation has sufficient freedom to evade this problem, 
provided that the DM scalar is specifically a real singlet. By choosing suitable values for the 
non-minimal couplings at the Planck scale, it is possible for the nnitarity-violation scale to be 
much larger than the inflaton field throughout inflation, so that the predictions of the model 
are robust. Therefore, in addition to providing a minimal candidate for WIMP DM, the 
extension of the SM by a non-minimally conpled real gauge singlet scalar can also account 
for inflation while having a consistent scale of nnitarity-violation. 

The paper is organized as follows. In section 2 we review the real gange singlet scalar ex- 

^The case of singlet DM added to Higgs Inflation was considered in [10]. 
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tension of the SM and the S-inflation model. We estimate the predictions of the model for the 
spectral index Hs and discuss the effect of constraints from inflation on the model parameter 
space. Section 3 considers the DM phenomenology of the model and the implications from 
DM searches. In section 4 we discuss how to connect these two aspects via renormalisation 
group evolution and which constraints follow from electroweak vacuum stability and pertur- 
bativity. The scale of unitarity-violation during inflation and the consistency of S-inflation 
are discussed in section 5. Finally, we present our results in section 6 and our conclusions in 
section 7. Additional details are provided in the Appendix. 


2 The S-inflation model 


S-inflation is a version of the non-minimally coupled inflation model of [7] in which the scalar 
held is identihed with the gauge singlet scalar responsible for thermal relic cold DM. In the 
present work, we focus on the case of a real singlet scalar s. In the Jordan frame, which is 
the standard frame for interpreting measurements and calculating radiative corrections, the 
action for this model is 


Sj 


^d^x\c^+{d^H)^ + -d^sd^^s 


V{s^, H^H) 


( 2 . 1 ) 


where is the SM Lagrangian density minus the purely Higgs doublet terms, mp is the 
reduced Planck mass and 


Vis^H^H) = Xh 


(h^h) 


+ 2 S ' 


'H^H+-Xs 
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( 2 . 2 ) 


with V = 246GeV the vacuum expectation value of the Higgs held. Writing H = {h+v, 0)/\/2 
with a real scalar h we obtain 

V{s‘^,h) = V (h) + ^ Xs ^ Xhs V h ^ Xhs , (2.3) 


where we have introduced the physical singlet mass + Xhs n^/ 2 . 

In order to calculate the observables predicted by inhation, we perform a conformal 
transformation to the Einstein frame, where the non-minimal coupling to gravity disappears. 
In the case that s 7 ^ 0 and h = 0, this transformation is dehned by 


Shi' 



= 1 + 


mt 


The transformation yields 
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where R is the Ricci scalar with respect to 5 ^,^. We can then rescale the held using 

dxs _ + s^/ml 

ds 94 
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which gives 


5e 



'^ + \9^'^d,Xsd,Xs-U{xs,^^ 


(2.7) 


with 

= (2,8) 

The relationship between s and Xs is determined by the solution to eq. (2.6). In particular, 
for s ^ mp/v^, the Einstein frame potential is 


C/(Xs,0) 


As rup 


1 + exp 


‘^Xs 

\/6mp 


-2 


(2.9) 


This is sufficiently flat at large to support slow-roll inflation. 

An analogous expression is obtained for the potential along the /i-direction. In both 
cases, the Einstein frame potential is proportional to where (j) = s or h. Therefore 

the minimum of the potential at large s and h will be very close to /i = 0 and inflation will 
naturally occur along the s-direction if <C Xh/^f^, which is true for example if 

and As ~ Aft. 

In the following, inflation is always considered to be in the direction of s with h = 0. 
The conventional analysis of inflation can then be performed in the Einstein frame. After 
inflation, the Jordan and Einstein frames will be indistinguishable since ^g <C rrip and 
so n —)• 1. Therefore the curvature perturbation spectrum calculated in the Einstein frame 
becomes equal to that observed in the physical Jordan frame at late times. 

The classical (tree-level) predictions for the spectral index and tensor-to-scalar ratio 
are [33] 


tree 

s 




0.965 , 


^tree 


^ + o 

N2 



= 3.6 X 10 “^ , 


( 2 . 10 ) 

( 2 . 11 ) 


while the field during inflation is 


~ 4m|iV/(36) ■ 


( 2 . 12 ) 


In the equations above N is the number of e-foldings as defined in the Einstein frame, 
which differs from that in the Jordan frame by iV ~ iV -|- ln(l/\/]V) [8], and we have used 
N = 58. ^ The classical predictions are in good agreement with the most recent Planck 
values, Ug = 0.9677 ± 0.0060 (68% confidence level (CL), Planck TT -|- lowP -|- lensing) and 
^ 0.002 < 0.11 (95% CL, Planck TT -|- lowP -|- lensing) [1]. 

The classical predictions for S-inflation are the same as those of any model based on 
the SBB model. Differences between S-inflation and other models do however arise from 

^Reheating in S-inflation occurs via stochastic resonance to Higgs bosons through the coupling Xhs- It 
was shown in [16] that this process is very efficient and makes quite precise predictions for the reheating 
temperature and the number of e-foldings of inflation, with 57 < iV < 60 at the WMAP pivot scale. This in 
turn allows for quite precise predictions of the inflation observables. 
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quantum corrections to the effective potential. To include these corrections, we calculate 
the RG evolution of the various couplings as a function of the renormalisation scale (see 
section 4 for details). We can then obtain the renormalisation group (RG)-improved effective 
potential for s in the Jordan frame by replacing the couplings in eq. (2.3) by the running 
couplings and setting fi equal to the value of the field. For h = 0 (and neglecting the singlet 
mass term, nig ~ 1 TeV <C s), this approach yields 

V^o{s\0) = . (2,13) 


The RG-improved potential can then be transformed into the Einstein frame in order to 
calculate the observables predicted by inflation.^ The inflationary parameters are calculated 
using the methods discussed in [16]. In particular, the Einstein frame slow-roll parameters 
are given by 


nip / 1 dU \ ^ 
nip d'^U 

:2 _ dU d^U 

1/2 • 


(2.14) 


3 Singlet scalar as dark matter 

Let us now turn to the phenomenology of the singlet scalar in the present Universe and 
at energies well below the scale of inflation [19, 24, 35]. Most importantly, the assumed 
Z 2 symmetry ensures the stability of the scalar, so that it can potentially account for the 
observed abundance of DM [2, 3]. If the mass of the singlet is comparable to the electroweak 
scale, the singlet is a typical WIMP, which obtains its relic abundance from thermal freeze- 
out. Indeed, at low energies, where the effects of the non-minimal coupling to gravity are 
negligible, our model becomes identical to what is often referred to as Higgs Portal Dark 
Matter [21, 22, 36], because all interactions of the singlet with SM particles are mediated 
by the Higgs. In this section we review the constraints on these models and determine the 
parameter space allowed by the most recent experimental results. In the process, we point 
out several discrepancies in the literature and resolve the resulting confusion. 

3.1 Relic abundance 

The calculation of the relic abundance of singlet scalars is discussed in detail in [19]. Three 
kinds of processes are relevant for the annihilation of singlets into SM states: annihilation 
into SM fermions, annihilation into SM gauge bosons and annihilation into two Higgs par¬ 
ticles. The first kind dominates as long as rUg < f^w-, while for larger masses the second 
kind gives the largest contribution. Notably, all of these processes can proceed via an s- 
channel Higgs boson, leading to a resonant enhancement of the annihilation cross-section 
and a corresponding suppression of the DM relic abundance for rug ~ mh(2.‘^ 

^In [34] it was proposed to use of the Einstein frame for the computation of quantum corrections. The 
Jordan frame analysis is, however, easier to implement correctly, being a straightforward extension of the 
Standard Model analysis. 

^The process ss — > hh also receives a contribution from t-channel singlet exchange, which gives a relevant 
contribution if \hs is large compared to Xh- 
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For the present work we calculate the singlet abundance using micrOMEGAs_3 [37], 
which numerically solves the Boltzmann equation while calculating the Higgs width in a 
self-consistent way. It is then straightforward to numerically find the coupling \hs that gives 
/i^ = 0.1197, in order to reproduce the value of the DM density /i^ = 0.1197 ± 0.0022 
determined by Planck (TT -|- lowP, 68% CL) [38]. For example, we find Xhs ~ 0.08 for 
nis = 300 GeV and X^s ~ 0.30 for = 1000 GeV. These values agree with the ones found 
in [19, 23], but disagree with [20-22, 24] by a factor of 2 after accounting for the different 
conventions.^ 


3.2 Direct detection constraints 

The strongest constraints on Xhs stem from DM direct detection experiments, since the 
singlet-Higgs coupling induces spin-independent interactions between the DM particle and 
nuclei. The scattering cross-section at zero momentum transfer is given by [19] 

*^31 = —i-4 2 ’ 

47r rrig 

where rUn is the neutron mass, //r = {'ms'mn)/{pis + rUn) is the reduced mass and /at is 
the effective Higgs-nucleon coupling.® In terms of the light-quark matrix elements the 
effective coupling can be written as 


(3.1) 



q=u,d,s 


(3.2) 


The values of can either be determined phenomenologically from baryon masses and 
meson-baryon scattering data or computed within lattice QCD. A comparison of the different 
methods was recently performed in [19] and we adopt their result of /tv = 0.30 for the effective 
coupling. 

The scattering cross-section given above can be directly compared to the bound obtained 
from the LUX experiment [18]. Indeed, as shown in figure 1, LUX is typically sensitive to the 
same range of values for Xhs as what is implied by the relic density constraint. Specifically, 
the LUX bound excludes the mass ranges 5.7 GeV < m* < 52.6 GeV and 64.5 GeV < rus < 
92.8 GeV. 


3.3 Invisible Higgs decays 

Direct detection experiments cannot constrain singlet scalars with a mass of a few GeV or 
less, since such particles would deposit too little energy in the detector to be observable. This 
parameter region can however be efficiently constrained by considering how the Higgs-singlet 
coupling Xhs would modify the branching ratios of the SM Higgs boson. The partial decay 
width for /i —>• ss is given by' 


r{h —)• ss) 


327r ruh V 


(3.3) 


®The Higgs-singlet coupling is called Xhs in [19], Xhs in [20], Xhss in [21, 22], Adm in [23] and 02 in [24]. 
The respective conventions are captured by Xhs = Xhs = Xhs/'^ = Xhss/2 = Adm/2 = 2 02 . 

®Note that [24] uses an approximate expression valid for rris S> rn„, such that fj,^ ~ rn^. 

^This equation agrees with [4, 19-22, 35] but disagrees with [23]. 
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Figure 1. Excluded parameter regions from LUX (red, dotted) and searches for invisible Higgs decays 
(blue, dashed) compared to the coupling implied by the relic density constraint (green, solid). 



This theoretical prediction can be compared to the experimental bound on invisible 
Higgs decays from the LHC. Direct searches for invisible Higgs decays in the vector boson 
fusion channel give BR(/i —)• inv) < 0.29 [39]. A somewhat stronger bound can be obtained 
from the observation that in our model there are no additional contributions to the Higgs 
production cross-section and no modifications of the partial decay widths of the Higgs boson 
into SM final states. Therefore the presence of an invisible decay channel leads to an overall 
reduction of the signal strength in visible channels. A global fit of all observed decay channels 
(combined with the bounds on invisible Higgs decays) then gives BR(/i —)• inv) <0.26 [40]. 

Crucially, the bound from invisible Higgs decays becomes independent of the singlet 
mass for rUg << mhj2. Invisible Higgs decays will therefore provide the strongest constraints 
for small singlet masses. Indeed, this constraint rules out the entire mass region where direct 
detection experiments lose sensitivity (see figure 1). Asa result, only two mass regions remain 
viable: a low-mass region 52.6 GeV < mg < 64.5 GeV and a high-mass region rrig > 93 GeV. 


3.4 Other constraints 

It has been pointed out recently [24] that bounds on y-ray lines from Fermi-LAT [41] rule 
out the parameter region where rrig is slightly above mhl2. To be safe from this constraint, 
we will focus on the mass range 52.6 GeV < rUg < 62.4 GeV, which we shall refer to as the 
low-mass region. For the high-mass region, on the other hand, there are no strong constraints 
from indirect detection. Moreover, collider searches for singlet scalars with rrig > mh/2 are 
extremely challenging [20-23, 35] and consequently, there are no relevant bounds from the 
LHC for the high-mass region [42]. 

The most significant improvements in sensitivity in the near future are expected to 
come from direct detection experiments. Indeed, XENONIT [43] is expected to improve 
upon current LUX constraints on the DM scattering cross section by a factor of about 50 
and will therefore be able to probe the high-mass region up to nrig ~ 4 TeV. As we will show, 
in S-inflation singlet masses larger than about 2 TeV are excluded by the Planck 2-a upper 
bound on rig and perturbativity. XENONIT will therefore be able to probe the entire high- 
mass region relevant for singlet inflation. Similarly, XENONIT can also further constrain the 
low-mass region and potentially probe singlet masses in the range 53 GeV < mg < 57 GeV. 
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4 Renormalisation gronp evolution and theoretical constraints 


In order to connect the inflationary observables of our model to the measured SM parameters 
and the DM phenomenology discussed in the previous section, we need to calculate the 
evolution of all couplings under the RG equations [8, 10, 12, 14, 44-47]. Existing analyses 
have considered the RG equations for the SM at two-loop order and the contributions of the 
singlet sector and non-minimal coupling at one-loop order (see also the Appendix) [8, 10]. 
To examine the issue of vacuum stability, we improve the accuracy of our analysis further by 
incorporating the three-loop RG equations for the SM gauge couplings [48] and the leading 
order three-loop corrections to the RG equations for and yt [49].^ 

When considering large field values for either s or h, the RG equations are modified, 
because there is a suppression of scalar propagators. This suppression is captured by inserting 
a factor 




1 + 


^ 4 > 

mp 




(4.1) 


with (f> = s {4> = h) for each s {h) propagating in a loop [8]. The changes in the RG equations 
for large values of the Higgs held have been discussed in detail in [47]. The modifications 
resulting from large singlet field values can be found in [8] and are reviewed in the Appendix. 
Note that, when considering S-inflation, such that s S> /i, we can set the suppression factor 
Ch = 1- 

We determine the values of the SM parameters at y = mt following [46]. Using the 
most recent values from the Particle Data Group [50] 


mt = (173.2 ± 0.9) GeV, = (125.09 ± 0.24) GeV, 05 ( 771 ^) = 0.1185 ± 0.0006 , (4.2) 


we obtain at y = nit 

yi = 0.936 ± 0.005, A/, = 0.1260 ± 0.0014, ys = 1.164 ± 0.003 . 


(4.3) 


Unless explicitly stated otherwise, we will use the central values for all calculations below. 
For given couplings at y = mt, we then use the public code RGErun 2.0.7 [51] to calculate 
the couplings at higher scales. 

In contrast to the remaining couplings, we fix the non-minimal couplings and at 
y, = nip. In order to obtain the correct amplitude of the scalar power spectrum, we require 

y = (0.00271 TTip)"^ , (4.4) 

where U and e are the potential and the first slow-roll parameter in the Einstein frame 
at the beginning of inflation, as defined in eq. (2.8) and eq. (2.14) respectively. Imposing 
equation (4.4) allows us to determine at the scale of inflation once all other parameters 
have been fixed. Note that, since the value of s at the beginning of inflation also depends on 
^s, equation (4.4) can only be solved numerically. We then iteratively determine the required 
value of at the electroweak scale such that RG evolution yields the desired value at the 
scale of inflation. 

The coupling plays a very limited role for the phenomenology of our model because we 
do not consider the case of large Higgs held values for inflation. As a result our predictions 

®We thank Kyle Allison for sharing his numerical implementation of these equations. 
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for the inflationary observables show only a very mild dependence on so that can 
essentially be chosen arbitrarily. Nevertheless, it is not possible to simply set this parameter 
to zero, since radiative corrections induce a mixing between and ^h- Moreover, we will 
see below that the value of plays an important role for determining whether our model 
violates unitarity below the scale of inflation. As with we fix at mp and then determine 
iteratively the required value of at the electroweak scale.® 

4.1 Metastability 

We now discuss various theoretical constraints related to the RG evolution of the parame¬ 
ters in our model. It is a well-known fact that for the central values of the measured SM 
parameters, the electroweak vacuum becomes metastable at high scales, because the quartic 
Higgs coupling Xh runs to negative values (see e.g. [46]). This metastability is not in any 
obvious way a problem, as the lifetime of the electroweak vacuum is well above the age of 
the Universe [52] (note, however, that this estimate may potentially be spoiled by effects 
from Planck-scale higher-dimensional operators [53-55]). However, one may speculate that 
a stable electroweak vacuum is necessary for a consistent theory, for example if the vacuum 
energy relative to the absolute minimum is a physical energy density leading to inflation. It 
is therefore an interesting aspect of singlet extensions of the SM that scalar singlets give a 
positive contribution to the running of Xh [56, 57]: 

+ 227^2 

In fact, it was shown in [58] that, for the case of a minimally-coupled singlet, Xhs can be 
chosen such that the electroweak vacuum remains stable all the way up to the Planck scale 
and at the same time (for appropriate choices of the singlet mass nig) the singlet obtains a 
thermal relic density compatible with the observed DM abundance.^® 

In order to study electroweak vacuum stability, we need to consider the potential in the 
/i-direction with s = 0.^^ Vacuum stability then requires that Xh{^j) > 0 for /r up to mp. 
In the present study, we consider both the case where Xhg is sufficiently large to stabilise 
the electroweak vacuum and the case where Xhs only increases the lifetime of the metastable 
vacuum, but does not render it completely stable. We focus throughout on the case where 
Xhs is positive. 

4.2 Examples 

Figure 2 shows an example for the evolution of scalar couplings (left) and the non-minimal 
couplings (right) under the RG equations discussed above. Solid lines correspond to the case 
s S> /i, which is relevant for inflation, while dotted lines correspond to /i ^ s, which is relevant 
for vacuum stability. We fix the scalar couplings at the weak scale, choosing Xhs = 0.25 and 
As = 0.01, such that the observed relic abundance can be reproduced for nig = 835 GeV. We 
consider As <C Xhg, in which case the value of .^s necessary to obtain the correct amplitude of 
the scalar power spectrum is reduced. For our choice, we find ^s ~ 10^. Note, however, that 

®The running of between mp and the scale of inflation is completely negligible, since the relevant 
diagrams are strongly suppressed for large values of s. 

^°Note that if the singlet mixes with the Higgs, there will be additional threshold effects at /r = rris from 
integrating out the singlet [59]. In the setup we consider, however, this effect is not important [60]. 

more detailed study of the potential in general directions with both s yf 0 and h yf 0 (along the lines 
of [61]) is beyond the scope of the present work. 
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A^^{mt)= 0 . 25 , As(m,)=0.01, ^h(mp)=1000 



Afisimt)=0.25, As{mt)=0.0^, ^i,(mp)=1000 



Figure 2. Running of the scalar couplings Xg and X^s (left) and of the non-minimal couplings 
and (right) as a function of the renormalisation scale fj, for a typical parameter point in the 
high-mass region. Solid lines show the running in the s-direction, while dotted lines correspond to the 
running in the h-direction. In the right panel, we also show the suppression factor Cg, which modifies 
the running in the s-direction at large field value. 


while Xfis exhibits only moderate running, Ag grows significantly with the renormalisation 
scale because its /3-function contains a term proportional to Choosing even smaller 
values of Ag at the weak scale will therefore not significantly reduce its value at the scale of 
inflation nor the corresponding value of 

An important observation from figure 2 is that X^ does not run negative and hence 
the electroweak vacuum remains stable all the way up to the Planck scale. The additional 
contribution from the singlet scalar is sufficient to ensure Xh > 10“^ for all renormalisation 
scales up to mp. For field values s > 10^® GeV S> h, the singlet propagator is suppressed, 
leading to a visible kink in the running of Xh- We show the propagator suppression factor 
Cg in the right panel of figure 2. One can clearly see how this suppression factor affects the 
running of which becomes nearly constant for s > 10^^ GeV. In this particular example, 
we have chosen ^/i(mp) = 1000 (for s ^ h). This choice, together with A* <C 1, implies that 
the running of ^s(/r) from rrit to mp is negligible. 

In the low-mass region, we are interested in much smaller values of Xhg, typically below 
10“^. A particular example is shown in figure 3 (left) for the representative choice Xhg{mt) = 
0.002 and Xg = 0.0005, which yields the observed relic abundance for rUg ~ 57 GeV. We 
observe that if Ag and Xhg are both small at the electroweak scale, these couplings exhibit 
only very little running up to the scale of inflation. For the same reason, it is impossible to 
influence the running of the Higgs couplings sufficiently to prevent Xhg from running negative 
at around 10^^ GeV. 

If Xhg is small, we can obtain the correct scalar power spectrum amplitude with a much 
smaller value of ^g during inflation. For the specific case considered in figure 3, we find 
^g ~ 10^. For these values of ^g and Xhg, the loop-induced corrections to ^h are very small 
and hence this coupling changes only very slightly under RG evolution. 

4.3 Perturbativity 

In order for our calculation of the running couplings and the radiative corrections to the 
potential to be reliable, we must require that all couplings remain perturbative up to the 
scale of inflation, which is typically I0^^~10^® GeV. This requirement is easily satisfied for 
the SM couplings, but needs to be checked explicitly for the couplings of the singlet, which 
can grow significantly with increasing renormalisation scale /r. We follow [58, 62] and use 
the requirement of perturbative unitarity to impose an upper bound on the scalar couplings. 
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Ahs(m()=0.002, As(mf)=0.0005, 



/j [GeV] 


IJ [GeV] 


Figure 3. Running of the scalar couplings Xh, and X^s (left) and of the non-minimal couplings 
and (right) as a function of the renormalisation scale /it for a typical parameter point in the low- 
mass region. In the right panel, we also show the suppression factor Cg, which modifies the running 
at large field value. Note that Xh runs negative for fi > 10^^ GeV. 


This procedure gives 

dvr 

Xs < — and Xhs < Stt . (4.6) 

As we will see below, the non-minimal coupling can be much larger than unity without 
invalidating a perturbative calculation. Nevertheless, if and are both very large, pro¬ 
cesses involving both couplings may violate perturbative unitarity, implying that there may 
be new physics or strong coupling below the scale of inflation. We will now discuss this issue 
in more detail. 

5 Unitarity-violation during inflation 

In this section we estimate the scale of perturbative unitarity-violation as a function of the 
background inflaton field. Note that by “unitarity-violation scale” we mean the scale at which 
perturbation theory in scalar scattering breaks down, so this may in fact indicate the onset 
of unitarity-conserving scattering in a strongly-coupled regime [28, 31]. We will consider 
unitarity-violation in the scattering of scalar particles corresponding to perturbations about 
the background field. In the case of a real scalar s and the fields of the Higgs doublet H, 
there are two distinct scattering processes we need to consider: (i) 6s hi -H- 6s hi and (ii) 
/ii/i 2 ■<->• /ii/i 2 , where 6s is the perturbation about the background s field and hi and /i 2 
are two of the Higgs doublet scalars. Scattering with the other scalars in H is equivalent 
to these two processes. We will use dimensional analysis to estimate the scale of unitarity- 
violation by determining the leading-order processes in the Einstein frame which result in 
unitarity-violating scattering. 

It will be sufficient to consider the Einstein frame Lagrangian for two real scalar fields 
ipi, where — using the notation of [30] — stands for either s or a component of the Higgs 
doublet. Unitarity-violation requires that there are two different scalars in the scattering 
process, since in the case of a single scalar there is a cancellation between s-, t- and u-channel 
amplitudes [63]. 

Since unitarity-violating scattering in the Jordan frame is due to graviton exchange via 
the non-minimal coupling to R, we cat set U = 0. The Einstein frame action for two real 
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scalars is then of the form 


5e = y cf^xv^ 


i<j 


ij 2 ^ 


(5.1) 


where 


Cii = - 




dn(l)i du4>i and Cij = 


6 ij (t>i 4>j 9^'" d^(f)i du(j)j 


m 


2 


(5.2) 


with 

02 = 1 + M. (5.3) 

mf, 

The interaction terms proportional to it ij are responsible for the dominant unitarity-violation 
in scattering cross-sections calculated in the Einstein frame. These interactions are the Ein¬ 
stein frame analogue of scalar scattering via graviton exchange in the Jordan frame due to 
the non-minimal coupling. To obtain the scale of unitarity-violation in terms of the physical 
energy defined in the Jordan frame, we first canonically normalize the fields in the Einstein 
frame, then estimate the magnitude of the scattering matrix element and finally transform 
the unitarity-violation scale in the Einstein frame back to that in the Jordan frame. 

In the following we will denote the inflaton by (pi (= s), which we expand about the 
background field, i.e. (pi = (pi + 6(pi. The Higgs doublet scalars are denoted by (p 2 (= hi) 
and </>3 (= /i 2 ). The corresponding canonically normalized scattering fields in the Einstein 
frame are then defined to be (^i, (p 2 and Once we have determined the interactions of 
the canonically normalised fields, we use dimensional analysis to estimate the scale of tree- 
level unitarity-violation. For this purpose we introduce appropriate factors of E to make 
the coefficient of the interaction terms in C dimensionless. Energy scales E in the Einstein 
frame are related to the ones in the Jordan frame via E = E/Q., where during inflation 

~ ii(pi/mp ^ N ^ 1. Unitarity conservation implies that the matrix element for any 
2 •(-)> 2 scattering process should be smaller than 0(1), so we can determine the scale of 
unitarity-violation (denoted by A) by determining the value of E that saturates this bound. 
This was demonstrated explicitly in [8], by comparing the dimensional estimate with the 
exact value from the full scattering amplitude. 

In the following, we consider three regimes for (pi, each of which leads to a different 
form of the Lagrangian and the scattering amplitudes: 

_2 

• Regime A: Large field values. In this regime we have H > 1 (implying that (pi/rrip > 
1) and 6if(pl/mp > 1. 

• Regime B: Medium field values. In this regime we have approximately H ~ 1 (implying 
that (p\/mp < 1), but still ^ii(t?i/mp > 1. 

• Regime C: Small field values. Finally we consider 6if (pi/mf, < 1, which in particular 
implies H ~ 1. 
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5.1 Regime A: Large field values 

In this case the canonically normalized fields are cpi = \/6mp and (^ 92,3 = The 

interaction leading to unitarity-violation in bs hi scattering is 

C D ^2^^4 (^1 + Sfpl) 4>2 d^Scpi du(j)2 . (5.4) 

77T/p 1 L 


This results in a 3-point and a 4-point interaction. After rescaling to canonically normalized 
fields, the 3-point interaction is 


I D 


mp 


</52 diJibipi du‘P2 ■ 


(5.5) 


This interaction can mediate </ 9 iV 92 ‘^W 2 scattering at energy E via ip 2 exchange, with a 

matrix element given dimensionally by |A4| ~ Unitarity is violated once |AI| ~ 1, 

therefore the unitarity-violation scale in the Einstein frame is 

Ag> ~ ^ . (5.6) 

where the superscript (3) denotes unitarity-violation due to the 3-point interaction. In the 
Jordan frame A) 2 ^ = UA) 2 ^ where 14 ~ (pilmp, therefore 

aS ~ • (5.7) 

Similarly, the 4-point interaction has an Einstein frame matrix element given by |AI| ~ 
^2 E^/rup, therefore the scale of unitarity-violation in the Jordan frame is 


A 


(4) 

12 


r\j 



(5.8) 


(3) (3) 

For ^2 > 1; this is larger than A) 2 ^ therefore A) 2 ' is the dominant scale of unitarity-violation. 
In general these estimates of the unitarity-violation scales are valid if the scalars can be 
considered massless, which will be true if 4>i < A) 2 ^ i.e. for > ^ 2 - 

In the case of Higgs scattering (/92<^3 ^ <^ 2923 , there is only the 4-point interaction 
following from 

C D ^ 2 ^^^ 02 03 g^'^ 9^02 0i.03 • (5.9) 

TTZp a L 

As the canonically normalized Higgs fields are in general given by (^ 2,3 = 02,3/^4, the unitarity- 
violation scale in the Einstein frame is generally A 23 = mp /^2 (since 6 = 6 T 

both scalars are part of the Higgs doublet). On translating the energy to the Jordan frame, 
the unitarity-violation scale becomes 

A23 ~ ^01 , (5.10) 

42 

(3) 

which is the same expression as for A) 2 b 
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5.2 Regime B: Medium field values 


In this case the canonically normalized fields are (pi = and (^2 = 4’2- Since 

n = 1, the energies are the same in the Einstein and Jordan frames. Using the same procedure 
as before, we find 


A 


(3) 

12 


r\j 


mp 

6 


(5.11) 


and 



(5.12) 


For Higgs scattering we obtain 


A23 ~ 


mp 

6 


(5.13) 


5.3 Regime C: Small field values 

In this case the Einstein and Jordan frames are completely equivalent. Therefore 


A 


(3) 

12 


mt 


66 ? 2</>1 


(5.14) 


and 


For Higgs scattering we obtain 


A 


(4) 

12 




mp 


A23 ~ 


mp 

6 


(5.15) 

(5.16) 


5.4 Discussion 

In summary, for > ih the smallest (and so dominant) scale of unitarity-violation in each 
regime is given by: 


A : A® 

B: 


Alt’ ~ 


mp 

/ 


(5.17) 


where A^^^ = A^ 2 ^ = ^ 23 , 6 = 6 and = 6- 

In the case of Higgs Inflation, the scale of unitarity-violation is obtained as above but 
with set equal to ^h, since the inflaton is now a component of H. During inflation A ~ 
with 4 r-u 10^. Since this energy scale is less than cpi, the gauge bosons become 
massive and only the physical Higgs scalar takes part in scattering. Since unitarity-violation 
requires more than one massless non-minimally coupled scalar, there is no unitarity-violation 
at energies less than (pi. Unitarity-violation therefore occurs at A ~ ss i.e. the 

unitarity-violation scale is essentially equal to the Higgs field value during inflation [14, 32]. 
As a result, either the new physics associated with unitarising the theory or strong coupling 
effects are expected to significantly modify the effective potential during inflation. It is 
uncertain in this case whether inflation is even possible and its predictions are unclear. 
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Af,simt)=0.25, As(/n/)=0.01, f/,(mp)=1000 


Ahs(m,)=0.25, As(m,)=0.01, ^h(mp)='\ 




Figure 4. The scale of unitarity-violation A as a function of the field value s (blue) for two different 
choices of parameters. The orange dashed line indicates the condition A > s, which must be satisfied 
in order to avoid unitarity-violation. In the left panel (with A/i(mp) = 500), unitarity is violated for 
s > 10^6 GeV. In the right panel (with Xh{mp) = 50), no unitarity-violation occurs up to the scale of 
inflation. Note that in the right panel ^/i(/i) runs negative for /i < 5 x 10^^ GeV. 


In S-inflation, on the other hand, it is possible to ensure that A S> s provided is 
sufficiently small compared to at the scale of inflation. This is illustrated in figure 4 for 
Chimp) = 1000 (left) and Chimp) = 1 (right), taking Xhsimt) = 0.25, Xsimt) = 0.01 and 
Csimp) ~ 10^ as above. Both panels show the scale of unitarity-violation A as a function 
of the field value s. For Chimp) = 1000 we observe that A < s for s > 10^® GeV, which is 
significantly smaller than the field value at the beginning of inflation. For Chimp) = 1, on 
the other hand, A always remains larger than s. In this case it is reasonable to assume that 
new physics, in the form of additional particles with mass of order A (or strong coupling 
effects^^), will have only a small effect on the effective potential at the scale /r = s. 

The right panel of figure 4 exhibits another new feature: We find A —)• oo for s ~ 
10^^ GeV. The reason is that, as already observed in figure 2, Chimp) exhibits a strong 
running for fi < 10^® GeV. Gonsequently, if we fix Ch to a rather small value at the Planck 
scale, e.g. Chimp) = 1, Chit^) will run negative at lower scales. During the transition, Chit^) 
will be very small and hence the scale of unitarity-violation can be very large. 

It should be emphasized that the advantage of S-inflation with respect to unitarity- 
violation is only obtained if the singlet is a real scalar. In the case of a complex singlet, the 
real and imaginary parts of the scalar both have the same non-minimal coupling Cs ■ Therefore 
we would have = C 2 in the above analysis and the unitarity-violation scale would become 
the same as in Higgs Inflation. Therefore the requirement that unitarity is not violated 
during inflation predicts that the DM scalar is a real singlet scalar. 

6 Results 

In this section we combine the experimental and theoretical constraints discussed above and 
present the viable parameter space for our model. Out of the five free parameters (As, Xhs, 
ms, Cs and Ch) we can eliminate X^g (or mg) by requiring the model to yield the observed 
DM abundance (see figure 1) and Cs by imposing the correct amplitude of the scalar power 

unitarisation by strong coupling, A is automatically field dependent and equal to the scale at which 
the potential is expected to change. In unitarisation by new particles, on the other hand, A is only an upper 
bound on the masses of the new particles. Moreover the masses need not be field-dependent in order to 
unitarise the theory. Therefore strong coupling is more naturally compatible with the scale of inflation. 
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Figure 5. Predictions for inflation in the high-mass region (500 GeV < rus < 2500 GeV) as a function 
of rUs and Ag for fixed = 100. Shown are the deviations from the tree-level predictions 

^tree _ Q 955 (left) as well as the value of at the beginning of inflation (right). For each value of nis 
the coupling Xfig has been fixed by the relic density requirement, as shown on the top of each panel. 
The grey shaded region indicates the parameter region where couplings become non-perturbative 
below the scale of inflation and the purple shaded region indicates the parameter space where Xh 
runs negative below the scale of inflation, leading to a metastable vacuum. We furthermore show 
the parameter region excluded by the upper bound on rig from Planck at 95% CL (shaded in blue). 
The green line in the right panel indicates the value of ^g where unitarity is violated at the scale of 
inflation. 


spectrum. In the following we always ensure that these two basic requirements are satisfied, 
and then consider additional constraints in terms of the remaining parameters Xg, mg (or 
Xhg) and We begin with a detailed discussion of the high-mass region and then turn to 
the low-mass region. 

6.1 The high-mass region 

Let us for the moment fix ^hi'nT'p) = 100 and study how the predictions depend on Ag and mg 
(or, alternatively, X^g)- The left panel of figure 5 shows the predicted value of Ug compared 
to the tree-level estimate = 0.965. We find that in our model Ug is always slightly 
larger than the tree-level estimate, but the differences are typically An* < 0.01. Only for 
Xhg > 0.5, corresponding to mg > 2 TeV, do the differences grow so large that the model can 
be excluded by the Planck 2-a bound, Ug < 0.98. In the same parameter region we find the 
largest differences between the tree-level predictions of r and a (see section 2) and the value 
predicted in our model. However, we find these deviations to be negligibly small. In particular 
our model predicts r < 0.01 everywhere, i.e. the tensor-to-scalar ratio would be very difficult 
to observe in the near future^^. Figure 5 also shows the parameter region excluded by the 
requirements that all couplings remain perturbative up to the scale of inflation (shaded in 
grey). This constraint requires Xg < 0.3 for small values of Xhg and becomes more severe 
with increasing Xhg- 

Finally, we also show the parameter region where Xhg is too small to prevent Xh from 
running to negative values in the /i-direction (shaded in purple). While this is not fatal for 

^^Next generation CMB satellites, such as PIXIE [64] and LiteBIRD [65], plan to measure r to an accuracy 
of 5r < 0.001. This would be sufficient to detect the tensor-to-scalar ratio in our model. 
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the model (the electroweak vacuum remains metastable with a lifetime that is longer than the 
one predicted for the SM alone), this constraint may be physically significant depending upon 
the interpretation of the energy of the metastable state. We note, however, that the bound 
from metastability depends very sensitively on the assumed values of the SM parameters 
at the electroweak scale. Indeed, it is still possible within experimental uncertainties (at 
95% CL) that the electroweak vacuum is completely stable even in the absence of any new 
physics [46]. 

For the currently preferred values of the SM parameters, we find the interesting pa¬ 
rameter region to be 0.2 < Xhs < 0.6 corresponding roughly to 700 GeV < nis < 2 TeV. 
Very significantly, this entire range of masses and couplings can potentially be probed by 
XENONIT. 

To study the predictions of inflation — and in particular the scale of unitarity-violation — 
in more detail, we show in the right panel of figure 5 the value of (at the scale of inflation) 
required by the scalar power spectrum amplitude. We typically find values around 10^, al¬ 
though values as large as 10^ become necessary as A* comes close to the perturbative bound. 
Since we have fixed ^/i(mp) = 100 in this plot, such large values of imply that \^/^h > Ij 
which in turn means that the scale of unitarity-violation is larger than the field value s at 
the beginning of inflation. Conversely, if both A* and Xhs are small, can be significantly 
below lO'^, such that unitarity is violated below the scale of inflation. The parameters for 
which the scale of unitarity-violation is equal to the scale of inflation is indicated by a green 
line. 

Let us now turn to the dependence of our results on the value of .^h(mp). For this 
purpose, we fix A* = 0.01 and consider the effect of varying ^/i(mp) in the range 0 < ^/i(mp) < 
1000. We find that neither the constraint from Planck, nor the bounds from metastability and 
perturbativity depend strongly on ^/i(mp). Nevertheless, as discussed in section 5, C/i(mp) 
does play a crucial role for the scale of unitarity-violation. We therefore show in the left panel 
of figure 6 the scale of unitarity-violation at the beginning of inflation divided by the field 
value s at the beginning of inflation. This ratio is to be larger than unity in order to avoid 
unitarity-violation. As indicated by the green line, this requirement implies ^h(mp) < 150 in 
the parameter region of interest. 

As discussed in section 4, small values of ^hi'mp) imply that ^/x(/i) will run to negative 
values for ^ ^ mt- To conclude our discussion of the high-mass region we therefore show 
in the right panel of figure 6 the magnitude of as a function of mg and ^/i(mp). The 

thick black line indicates the transition between ^hii^^t) > 0 and ^^(mi) < 0. By comparing 
this plot with the one to the left, we conclude that within the parameter region that avoids 
unitarity-violation ^h(?^t) necessarily becomes negative. 

6.2 The low-mass region 

We study the predictions for the low-mass region in figure 7. The top row shows An^ and 
at the beginning of inflation as a function of A<j and Xhs for C/i(mp) = 1. These plots 
are analogous to the ones for the high-mass region in figure 5. The crucial observation 
is that, unless A^ <C Xhs, radiative corrections to the inflationary potential are completely 
negligible, because any contribution proportional to A^ is suppressed by powers of c* <C 1 
during inflation. Consequently, in most of the low-mass region, Ug is identical to its tree-level 
value. We furthermore find that, as expected, the tensor-to-scalar ratio and the running of 
the spectral index are both unobservably small. Since radiative corrections play such a small 
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Figure 6. Left: the scale of unitarity-violation A compared to the field value at the beginning of 
inflation s as a function of ^hijn^p) and Ws (or Xus) for As = 0.01. In order to avoid unitarity-violation, 
we require log A/s > 0. Right: the corresponding value of at the electroweak scale = mt). The 
shaded regions correspond to the same constraints as in figure 5. 


role, the value of required from inflation depends almost exclusively on A^. As a result, it 
is easily possible to have < 1000 during inflation for As < 10 ^ and ^s < 100 for As < 10 

Since ^s can be much smaller in the low-mass region than in the high-mass region, it is 
natural to also choose a very small value for In fact, for typical values in the low-mass 
region ^s ~ 10^ and As ~ Xhs ~ 10“^, the loop-induced contribution to is < 10“^, so 
that it is technically natural to have <C 1. One then obtains ^ 1 and hence the 

scale of unitarity-violation is well above the scale of inflation. It is therefore possible without 
difficulty to solve the issue of unitarity-violation in the low-mass region. 

This conclusion is illustrated in the bottom row of figure 7, which should be compared 
to figure 6 from the high-mass case, except that we keep = 1 fixed and vary As 

instead. The bottom-left plot clearly shows that (for our choice of ^h) the scale of unitarity- 
violation is always well above the field value at the beginning of inflation, so that the problem 
of unitarity-violation can easily be solved in the low-mass region. In addition, if the scalar 
couplings are sufficiently small, will have negligible running from the Planck scale down to 
the weak scale. It is therefore easily possible to set e.g. ^hi'mp) = 1 and still have > 0, 

as illustrated in the bottom-right panel of figure 7. It is not possible, however, to ensure at the 
same time that Xh remains positive for large field values of /i. In other words, the electroweak 
vacuum is always metastable in the low-mass region (for the preferred SM parameters). 

In the high-mass region we found that the coupling Xhs is bounded from below by the 
desire to stabilise the electroweak vacuum and from above by constraints from Planck and 
the requirement of perturbativity. In the low-mass region, on the other hand, we obtain an 
upper bound on Xhs from LUX and a lower bound on Xhs from the relic density requirement. 
Compared to the high-mass region, the allowed range of couplings in the low-mass region is 
much larger and therefore much harder to probe in direct detection experiments. If indeed 
nis is very close to mh/2, and As, Xhs < 10“^, it will be a great challenge to test the model 
predictions with cosmological or particle physics measurements. 
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Figure 7. Top row: Predictions for inflation in the low-mass region (52.5 GeV < mg < 62.5 GeV) as 
a function of the couplings \hs and As for fixed = 1. Shown are the deviations from the tree- 

level predictions ns"^®® = 0.965 (left) and the value of ^s at the beginning of inflation (right). Bottom 
row: The scale of unitarity-violation A compared to the held value at the beginning of inflation Sinf 
(left) and the value of at the electroweak scale (right). In all panels the grey shaded region indicates 
the parameter region where couplings become non-perturbative below the scale of inhation and the 
light blue shaded region represents the 95% GL bound from LUX. We do not show the metastability 
bound, since it covers the entire low-mass region. 


7 Conclusions 

The origin of inflation and the nature of DM are two of the fundamental questions of cos¬ 
mology. In the present work, we have revisited the possibility that both issues are unified 
by having a common explanation in terms of a real gauge singlet scalar, which is one of the 
simplest possible extensions of the SM. Considering the most recent experimental constraints 
for this model from direct detection experiments, the LHC and Planck, we have shown that 
large regions of parameter space remain viable. Furthermore, we find that in parts of the 
parameter space the scalar singlet can stabilise the electroweak vacuum all the way up to the 
Planck scale, while at the same time avoiding the problem of unitarity-violation present in 
conventional models of Higgs inflation. 

The scalar singlet can efficiently pair-annihilate into SM particles via the Higgs portal, 
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so that it is straight-forward in this model to reproduce the observed DM relic abundance 
via thermal freeze-out. We find two distinct mass regions where the model is consistent with 
experimental constraints from LUX, LHC searches for invisible Higgs decays and Fermi-LAT: 
the low-mass region, 53 GeV < < 62.4 GeV, where DM annihilation via Higgs exchange 

receives a resonant enhancement, and the high-mass region, > 93 GeV, where a large 
number of annihilation channels are allowed. 

In both mass regions it is possible without problems to fix the non-minimal couplings 
and in such a way that inflation proceeds in agreement with all present constraints. 
In particular, the tensor-to-scalar ratio and the running of the spectral index are expected 
to be unobservably small. On the other hand, radiative corrections to the spectral index 
typically lead to a value of Ug slightly larger than the classical estimate, i.e. Ug > 0.965. This 
effect is largest for large values of nig and Xhg and current Planck constraints already require 
nig < 2 TeV. The entire high-mass region compatible with Planck constraints will therefore 
be tested by XENONIT, which can constrain gauge singlet scalar DM up to nig ~ 4 TeV. 

In the high-mass region, the value of ^g required to obtain a sufficiently flat potential 
during inflation is typically ^g ~ 10^-10®. In spite of such a large non-minimal coupling, 
it is possible to have unitarity-conservation during inflation, in the sense that the scale of 
unitarity-violation can be much larger than the inflaton field. The reason is that only ^g 
needs to be large in order to reproduce the observed density perturbation, while the Higgs 
non-minimal coupling can be arbitrarily small. In the limit —?■ 0 there will be only one 
non-minimally coupled scalar field and therefore no unitarity-violation, provided that the 
inflaton, and so the DM particle, is a real scalar. 

We find that at large singlet field values s the scale of unitarity-violation is given by 
A ~ sy/^g/^h- If the non-minimal couplings satisfy .^s(mp) » ihi'm-p) at the Planck scale, it 
is possible for the unitarity-violation scale during inflation to be orders of magnitude larger 
than s. Such a hierarchy of couplings is stable under radiative corrections and consistent with 
the assumption that inflation proceeds along the s-direction. Furthermore, in the low-mass 
region X^g and can be so small that ^g ~ 10^-10^ is sufficient to obtain a flat enough 
potential. 

We conclude that it is possible for the inflaton potential in S-inflation to be safe from 
new physics or strong-coupling effects associated with the unitarity-violation scale. This 
contrasts with the case of Higgs Inflation, where unitarity is always violated at the scale of 
the inflaton. 

Another interesting observation is that if the singlet mass and the coupling X^g are 
sufficiently large (roughly nig > 1 TeV and Xhg > 0.3), the presence of the additional scalar 
singlet stabilises the electroweak vacuum, because the additional contribution to pre¬ 
vents the quartic Higgs coupling from running to negative values. This observation becomes 
important if a metastable electroweak vacuum is physically disfavoured, for example if the 
potential energy relative to the absolute minimum defines an observable vacuum energy. 

Given how tightly many models for DM are constrained by direct detection and LHG 
searches and how strong recent bounds on models for inflation have become, it is quite 
remarkable that one of the simplest models addressing both problems still has a large allowed 
parameter space. Nevertheless, the model is highly predictive. In particular, if the DM scalar 
is also the inflaton and unitarity is conserved during inflation, then DM is predicted to be a 
real scalar. Direct detection experiments will soon reach the sensitivity necessary to probe 
the entire parameter space relevant for phenomenology, with the exception of a small window 
in nig close to the Higgs resonance. The next few years will therefore likely tell us whether 
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indeed a singlet scalar extension of the SM can solve two of the central problems of particle 
physics and cosmology. 
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Appendix: RG Equations 

The RG equations for the scalar couplings can be obtained using the techniques detailed 
in [66-68], as in [8]. The one-loop /3-functions for the scalar couplings are 


= -^Vt + ^ ( 9 ^ + 5 '^)^) + -3g'^ 12 yj) Xt 

+ (I 8 c\ + 6) Xl + - cl Xlg , 


(A-1) 


16vr^ — AchCs -|- 6 (c| -|- l) Xh Xhs ~ 2 


+ 6 Xhs + 6 Cg As Xhs ^ 


(A-2) 


Ihvr^ + ^)^hs + 18 Cs As . 


(A-3) 


The propagator suppression factors are given by 



(A-4) 


where 4> is s or h. 


The RG equations for the non-minimal coupling can be derived following [69], as in [8] 
(see also [10]). One obtains 



(A-5) 
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